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ABSTRACT 

Aims. A huge multiwavelength campaign targeting the blazar AO 0235+164 was organized by the Whole Earth Blazar Telescope 
(WEBT) in 2003-2005 to study the variability properties of the source. 

Methods. Monitoring observations were carried out at cm and mm wavelengths, and in the near-IR and optical bands, while 
three pointings by the XMM-Newton satellite provided information on the X-ray and UV emission. 

Results. We present the data acquired during the second observing season, 2004-2005, by 27 radio-to-optical telescopes. The 
~ 2600 data points collected allow us to trace the low-energy behaviour of the source in detail, revealing an increased near-IR 
and optical activity with respect to the previous season. Increased variability is also found at the higher radio frequencies, down 
to ~ 15 GHz, but not at the lower ones. While the X-ray (and optical) light curves obtained during the XMM-Newton pointings 
reveal no significant short-term variability, the simultaneous intraday radio observations with the 100 m telescope at Effelsberg 
show flux-density changes at 10.5 GHz, which are more likely due to a combination of intrinsic and extrinsic processes. 
Conclusions. The radio (and optical) outburst predicted to peak around February-March 2004 on the basis of the previously 
observed 5-6 yr quasi-periodicity did not occur. The analysis of the optical light curves reveals now a longer characteristic time 
scale of variability of ~ 8 yr, which is also present in the radio data. The spectral energy distributions corresponding to the 
XMM-Newton observations performed during the WEBT campaign are compared with those pertaining to previous pointings of 
X-ray satellites. Bright, soft X-ray spectra can be described in terms of an extra component, which appears also when the source 
is faint through a hard UV spectrum and a curvature of the X-ray spectrum. Finally, there might be a correlation between the 
X-ray and optical bright states with a long time delay of about 5 yr, which would require a geometrical interpretation. 

Key words, galaxies: active - galaxies: BL Lacertae objects: general - galaxies: BL Lacertae objects: individual: AO 0235+164 
- galaxies: jets - galaxies: quasars: general 



1. Introduction 

Among the active galactic nuclei, blazars (i.e. flat- 
spectrum radio quasars and BL Lac objects) exhibit the 
most extreme properties: violent variability at all wave- 
Send offprint requests to: C. M. Raiteri, e-mail: 
raiteri@to . astro . it 
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lengths, high optical (and radio) polarization, superlu- 
minal motion of radio components, and brightness tem- 
peratures often exceeding the inverse Compton limit. 
These features are explained by assuming that the emis- 
sion comes from a relativistic plasma jet closely aligned 
with the line of sight (see e.g. iBlandford fc Reesl Il978t 
iBlandford fc Konigj flOTflh . 

One of the most intriguing blazars is AO 0235+164 
(J0238+1636) at z = 0.94 which, in addition to the above 
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characteristics, also shows a complex environment, where 
a couple of foreground galaxies at z = 0.524, possibly in- 
teracting, can both significantly contaminate its emission 
in the optical-UV band and absorb its radiation from the 
near-IR to the soft X-rays. Moreover, the stars of these 
objects might act as microlenses . producing at least part 
of the observed variability (e.g. lOstriker fc Vietri 1 19851 
IStickel et alllTflRa. 

The analysis of the radio and optical light curves of AO 
0235+164 from 1975 to 2000 revealed correlated variabil- 
ity and a quasi-periodicity of the ma in radio (and optical ) 
outbursts on a ~ 5.7 yr time scale ifRaiteri et all lioOll) . 
The latter feature made this obje ct a good candidate 
for the binary black hole scenario ^Romero et all [2003: 
lOstorero et all 12004 1 ), as in the cas e of OJ 287, another 
famo us blazar (e.g. ISillannaa et all Il988t IVillata et all 

In 2003-2005 AO 0235+164 has been the target of 
a huge multiwavelength observing effort chiefly aimed at 
verifying its possible quasi-periodicity, since a radio out- 
burst had been predicted to peak around February-March 
2004 (with an uncertainty of 6 months). The Whole Earth 
Blazar Telescope 1 (WEBT) consortium organized a radio, 
near-IR, and optical long-term monitoring campaign, dur- 
ing which three observations by the XMM-Newton satel- 
lite were performed, with simultaneous dense radio moni- 
toring at the 100 m Effelsberg telescope. A spectroscopic 
monitoring was also carried out at the VLT, TNG, and 
NTT telescopes to study the relationship between the 
Mgn line flux and the continuum flux density, and 15 
VLBA epochs were granted to study the source structure 
variability. The results of these observations will be pre- 
sented elsewhere. 

The results of the first obs e rving season 2003-2004 
were reported by iRaiteri et alJ l(2005h . while a detailed 
analysis of the X- ray observations was performed by 
IRaiteri et all (|2006[) . Here we present the radio-to-optical 
data acquired by 27 telescopes during the second observ- 
ing season, 2004-2005. The data of this extensive WEBT 
campaign 2003-2005, added to the historical light curves, 
allow one to study the behaviour of this source over the 
last 30 yr. 

This paper is organized as follows: the optical, near-IR, 
and radio light curves obtained during the 2004-2005 ob- 
serving season are presented in Sect. 2. In Sect. 3 we derive 
optical colour indices over a longer time interval, includ- 
ing the 1989-1991 and 1997-1998 outbursts, to study the 
spectral variability of the source. The intraday observa- 
tions performed with the Effelsberg radio telescope during 
the August 2004 and January 2005 XMM-Newton point- 
ings are presented and discussed in Sect. 4. The long-term 
optical and radio behaviour of AO 0235+164 is presented 
in Sect. 5. In Sect. 6 we analyse 11 SEDs of the source 
corresponding to different epochs of X-ray observations, 

1 http : / /ww w, to . a st ro . i t/blazars/webt/ see e.g. 
IVillata et all J2004al) ; iBottcher et alJ J2005I) : IRaiteri et all 
J2005I) . 



Table 1. Ground-based observatories participating in the 
2004-2005 observing season of the WEBT campaign. 



Observatory 


Tel. size 


Bands 


iVobs 


SAO RAS (RATAN) 


600 m a 


2.3, 4.8, 7.7, 


33 






11.1, 21.7 GHz 




Metsahovi 


14 m 


37 GHz 


142 


Noto 


32 m 


5, 8, 22 GHz 


56 


Medicina 


32 m 


5, 8 GHz 


23 


Effelsberg 


100 m 


4.9, 10.5 GHz 


68 


UMRAO 


26 m 


4.8, 8.0, 14.5 GHz 


112 


Pico Veleta 


30 m 


1, 2, 3 mm 


18 


Campo Imperatore 


110 cm 


J,H,K 


154 


Bohyunsan 


180 cm 


V,R,I 


126 


Mt. Maidanak 


150 cm 


U, B, V, R, I 


733 


Abastumani 


70 cm 


R 


328 


Crimean 


70 cm 


B,V,R,I 


100 


Skinakas 


130 cm 


B,V,R,I 


83 


Tuorla 


103 cm 


R 


20 


Armenzano 


40 cm 


B,V,R,I 


87 


Michael Adrian 


120 cm 


R 


48 


Asiago 


182 cm 


B,V,R,I 


38 


Torino 


105 cm 


R,I 


20 


Roque (KVA) 


35 cm 


R 


69 


Roque (NOT) 


256 cm 


U, B, V, R, I 


24 


Roque (TNG) 


358 cm 


U, B, V, R, I 


10 


CASLEO 


215 cm 


B,V,R,I 


108 


La Silla (NTT) 


358 cm 


R 


5 


Boltwood 


40 cm 


R 


129 


Bell 


60 cm 


R 


23 


Kitt Peak (RCT) 


130 cm 


R 


24 


Kitt Peak (WIYN) 


90 cm 


B,V,R,I 


12 


Total 






2593 


a Ring telescope 



trying to understand how many emission components are 
involved. Finally, the conclusions are drawn in Sect. 7. 

2. Observations during the 2004-2005 season 

The second part of the WEBT campaign took place 
from May 1, 2004 (JD = 2453127) to May 1, 2005 
(JD = 2453492). Table □ shows the list of the participat- 
ing observatories, with the size of the telescope, the fre- 
quency/wavelength/band, and the number of observations 
performed. As one can see, 2593 data points were acquired 
from 27 telescopes: 6 of them working at cm wavelengths, 
one at mm wavelengths, one in the near-IR bands, and 19 
in the optical bands. 

Since the source was often very faint, a number of data 
were affected by large errors. Hence, when constructing 
the light curves some data from the same telescope were 
binned 2 to increase the signal to noise ratio, while some 
other data, in disagreement with the trend of the light 
curves, were discarded. This cleaning process was mainly 
performed on the optical light curves, where the compari- 
son among the source behaviour in different bands and the 
better sampling can help to recognize unreliable points. 

2 With binning interval not exceeding half an hour. 
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Fig. 1. UBVRI light curves of AO 0235+164 during the observing season 2004-2005; vertical lines indicate the XMM- 
N ewton pointings of A ugust 2, 2004 and January 28, 2005. The ELISA contribution has been subtracted as explained 
in lRaiteri et all l)2005|) . 



Moreover, since the presence of an AG N only 2 arc- The UBVRI light curves of the 2004-2005 observing 
sec s outh of the source (named ELISA by iRaiteri et alJ season are shown in Fig. ^ Different symbols and colours 
2005) affects the optical photometry of AO 0235+164 es- are used to distinguish the provenance of the data, as de- 
pecially in the blue part of th e spectrum, we subt racted scribed in the legend. Vertical lines are drawn through the 
its contribution as discussed in lRaiteri et al.l 1 2005h . dates of the XMM-Newton pointings of August 2, 2004 

and January 28, 2005. 

With respect t o the previous observing season (see 
2.1. Optical and near-IR light curves IRaiteri et alJl200fil and the top panel of Fig. 0), in 2004- 

2005 the source showed increased variability, spanning 
Optical and near-IR data were collected as instrumental more than 2.5 mag, although the average level has re- 
magnitud es of the source and of the reference stars in mained rather faint. In the i?-band light curve, the best 
its field l|Smith et allll98"5llFiorucci et al.lll998|) in order sampled one, four major flares can be recognized (around 
to apply the same magnitude calibration to all datasets, JD = 2453210,2453250,2453310, and 2453390), which 
when possible. brought the source to be brighter than R = 17. The 
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Fig. 2. J UK light curves of AO 0235+164 during the ob- 
serving season 2004-2005 compared to the i?-band one 
(top). Vertical lines are drawn to guide the eye through 
variations in different bands. 



most dramatic variations were observed during the rising 
and dimming phases of the second flare, when the i?-band 
brightness first increased by ~ 1.0 mag in 2 days and then 
decreased by ~ 2.2 mag in 12 days. 

Near-IR monitoring was performed at the Campo 
Imperatore Observatory in the JHK bands. The resulting 
light curves are plotted in Fig. [2 compared to the i?-band 
one. Although the near-IR light curves are less sampled 
than the i?-band one, a general agreement is visible. In 
particular, it is possible to recognize the first flare around 
JD = 2453210 and the dip at JD = 2453350. 

2.2. Radio light curves 

A comparison between the i?-band light curve and the ra- 
dio behaviour at different wavelengths is displayed in Fig. 
H for both the observing seasons 2003-2004 and 2004- 
2005 of the WEBT campaign. Data from the VLA/VLBA 
Polarization Calibration Database are also shown 3 . The 
vertical lines mark the XMM-Newton pointings. Cubic 
splines through the 30-day binned light curve at 37 GHz 
and 50-day binned light curves at 22, 14.5, 8, and 5 



GHz are drawn to highlight the long-term variations 4 . 
Horizontal lines show the magnitude and flux-density av- 
erage values over the considered period: <R> = 17.84, 
<F 37 > = 1.66 Jy, <F 22 > = 1.32 Jy, <F 14 . 5 > = 1-44 Jy, 
<F 8 > = 1.42 Jy, <F 5 > = 1.49 Jy. These 2003-2005 op- 
tical and radio levels are relatively low if compared with 
those reached during the major historical outbursts; in- 
deed, when looking at the historical light curves (Fig. 
0) one can see that i? m i n = 14.03, F 3 7 jmax = 6.88 Jy, 
F 2 2,max = 6.10 Jy, Fi4. 5jmax = 8.09 Jy, F s , max = 7.03 Jy, 
and F 5 , max = 4.28 Jy. 

As already found for AO 0235+164 and for other 
blazars, variations have larger flux amplitudes at the 
higher frequencies. It is interesting to notice that a higher 
flux density is visible at 37 GHz at the beginning of the 
2004-2005 observing season, when the source was experi- 
encing an optical flare. The radio flux density from 37 to 
14.5 GHz was also higher at the time of the optical flare 
which occurred around JD = 2453300-2453330, while at 
the start of the 2003-2004 season some radio activity at 
all frequencies seems to follow the bright optical state. 
However, it is hard to say whether these higher radio 
states have some correlation with their contemporaneous 
optical events, or if they are delayed counterparts of pre- 
vious optical flares or, finally, if there is no correlation at 
all. Actually, each of these possibilities has already been 
singled out in the h istorical behaviour of the source (e.g. 
iRaiteri et all I2005L and references therein) . 



3. Optical colour indices 

The colour behaviour of blazars can be an important tool 
to disentangle different components contributing to the 
optical emiss i on va r iability . In the case of BL Lacertae, 
IVillata et all l|2002l l2004b|) showed that the short-term 
variations are strongly chromatic, with a bluer-when- 
brighter behaviour, likely due to intrinsic energetic pro- 
cesses. On the contrary, the longer-term flux changes do 
not imply important spectral variations, and are possi- 
bly of geometric origin. A recent stu dy of the quasar-type 
blazar 3C 454.3 l|Villata et all 12006^ shows that the opti- 
cal spectrum steepens when the source is in outburst and 
the jet contribution to the optical flux is dominant. During 
fainter optical states instead, the emission from a blue 
thermal disc may be responsible for the observed spec- 
tral flattening (bluer- when- fainter behaviour). The multi- 
colour optic al observa t ions o f 8 blazars (3 quasars and 5 
BL Lacs) bvlGu" et al.1 l|2006l) support this result, showing 
that quasars tend to be redder when they are brighter, 
while BL Lacs tend to be bluer. An analysis of the spec- 
tral behaviour of a s ample of blazars was performed by 
iFiorucci et alJ l|2004|k in particular, they found that the 
mean spectral slope is flatter during the rising phases of 
flares, while it is steeper during the declines. 



http : //www. vla.nrao . edu/astro/ calib/polar/ 



4 We chose a shorter time bin for the 37 GHz data because 
the light curve at 37 GHz is better sampled than those at the 
lower frequencies and flux density variations are faster. 
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Fig. 3. i?-band magnitudes (top) and radio fluxes (Jy) at different frequencies during the 2003-2005 WEBT campaign. 
Vertical lines mark the times of the three XMM-Newton pointings. Cubic splines on the radio data at 37, 22, 14.5, 8, 
and 5 GHz are plotted to better see the long-term trend. Horizontal lines indicate the mean fluxes at these frequencies. 
Two different symbols are used for the 37 GHz data from Metsahovi since they come from different datasets. 



We show the B - R colour index of AO 0235+164 as 
a function of time in Fig. 0] (bottom), compared to the R- 
band light curve (top). The plot covers the period 1989- 
2005 to include the last major optical outbursts. The B— R 
points have been obtained by coupling each B magnitude 
from a given telescope with the average of the R data 
from the same telescope acquired next the B datum (i.e. 
within 10-20 minutes). Only measurements with errors 
not greater than 0.08 mag were considered. In this way, 
we coupled 401 B data with 670 R ones to obtain indices 
with a maximum error of 0.10 mag. 

The colour index shows a noticeable scatter, from 1.27 
to 2.07, which implies important spectral variations. The 
average value is 1.73, with a standard deviation of 0.11. 
However, it is not possible to recognize any simple cor- 
relation with brightness. In particular, when considering 
the indices corresponding to the outburst phases of 1989- 
1991 and 1997-1998, they do not show peculiar values, 
since similar B — R values are found also in much fainter 
states. Thus, it seems that the long-term behaviour of the 
source is essentially achromatic. The V — I colour indices 
confirm this finding. 

However, one can notice a progressive spectral redden- 
ing in the last years: although the average R magnitude 5 in 
the 2002-2003, 2003-2004, and 2004-2005 observing sea- 
sons is similar (~ 17.8), the average B — R values are 1.66, 
1.72, and 1.82, respectively. 

Finally, one can derive the optical spectral index a op t 
(F„ oc v~ a ) from the colour index: a opt = \(B — R) — 
[A B - A R ) - 0.3521J/0.4125. By setting the extinction 
values Ab and Ar equal to 1 .904 and 1.260, respectively 
(see Table 5 in lRaiteri et al.ll2005l) . we obtain an average 
value a op t = 1.76, with standard deviation 0.26. Although 
this value has to be taken with some caution because it 
is calculated from two bands only, it nonetheless gives an 
indication of the optical spectral behaviour of the source. 
The spectrum appears rather steep, but still within the 
range of s pectral indices of low -energy peaked BL Lacs 
derived bv lFiorucci et al.l l)2004|) . The very high values of 
a op t (> 2) estimated for AO 0235+164 by various authors 
in previous works were obtained by neglecting the absorp- 
tion due to the z = 0.524 intervening system. 



4. Intraday radio observations at Effelsberg 

As in the case of the X MM-Newton pointin g of January 
18-19, 2004 discussed in lRaiteri et all l|2005l) . also during 
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the other two pointings of August 2, 2004 and January 
28, 2005, densely sampled radio observations were per- 
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Fig. 5. Radio light curves of the source and calibrators 
at 10.5 GHz (top, blue symbols) and 4.9 GHz (bottom, 
red symbols) obtained with the 100 m radio telescope at 
Effelsberg on August 2, 2004. The two calibrators are 3C 
067 (triangles) and 4C +09.11 (crosses). The horizontal 
lines indicate deviations of 1% from the mean flux den- 
sities, while the yellow areas highlight the period of the 
X-ray observations. 



formed at 4.9 and 10.5 GHz with the 100 m radio telescope 
at Effelsberg. The obse rving strategy and data reduction 
procedures are based on l Kraus et all l|2003f) and have been 
outlined in lRaiteri et al.l (|2005^l . The resulting light curves 
are shown in Fig. for the August 2004 pointing and in 
Fig.[5]for the January 2005 one. The yellow areas highlight 
the period of the X-ray observations 6 . Normalized flux 
densities of the two calibrators 3C 067 (J0224+2750) (tri- 
angles) and 4C +09.11 (J0238+1010) (crosses) are shown 
in the bottom panels. Horizontal lines in both the source 
and the calibrators panels indicate 1% deviations from 
the mean flux densities. They allow to grab the amount 
of variability and to compare the source one with that of 
the calibrators. 



5 Obtained by considering only the points which are used to 
calculate the colour indices. 



6 The X-ray light curves correspondin g to the three XMM- 
Newton pointings have been analysed in lRaiteri et alJ 12006): 
they show no significant variability. Furthe r analysis of the 
radio and X-ray light curves can be found in iKadleil J2005I) . 
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Fig. 4. B — R colour indices (bottom) compared to the i?-band light curve (top); in the top panel blue (larger) dots 
indicate the R data used to derive the colour indices. Horizontal lines are drawn through average values; vertical lines 
indicate the times of the X-ray satellite pointings, discussed in Sect. 5. 



In both epochs the 4.9 GHz light curve does not show 
significant variability, while the flux density changes at 
10.5 GHz are more interesting. Indeed, on August 2, 2004 
one can see first a sharp flux decrease of 3% in about 2 
hours 7 , and then a slower flux increase. On January 28, 
2005 the percent variation of the source flux is not greater 
than the calibrators one, but the fact that we see a quasi- 
linear trend suggests that this can be a piece of a larger 
variation, and that the probability that it is noise is very 
low. 

Table [21 reports some statistics: in both epochs the 
standard deviation a is greater than the root mean square 

uncertainty e = \J^2d = i ef ~~/N (where ej are the individ- 
ual errors) at 10.5 GHz, while it is smaller at 4.9 GHz. 
The same conclusion is reached by considering the reduced 
chi-square Xri which tests the hypothesis that the light 
curve can be modelled by a constant function: the values 
in Table|21tell that this hypothesis is fair for the 4.9 GHz 
data, but not for the 10.5 GHz ones. Finally, Table[2]also 
shows the modulation index m [%] = 100 a/<F u >, which 
provides a measure of the strength of the observed vari- 
ation amplitudes without taking into account the errors 
of the individual measurements. This is lower at 4.9 GHz 
than at 10.5 GHz and inferior to what was found for the 

7 Calculated as the variation between the first (the highest) 
and fifth (the lowest) point, which has also a slope equal to 
that determined by the second-fourth points. A similar, but 
more uncertain, slope is also simultaneously present in the 4.9 
GHz data. 



Table 2. Statistics on the Effclsberg data 



Freq. 


<F V > a e 


xi 


m 


(GHz) 


(Jy) (Jy) (Jy) 




(%) 




August 2, 2004 






10.5 


1.465 0.013 0.006 


3.76 


0.89 


4.9 


1.418 0.004 0.008 


0.29 


0.28 




January 28, 2005 






10.5 


1.250 0.012 0.009 


1.66 


0.96 


4.9 


1.311 0.005 0.007 


0.45 


0.38 



January 18-19, 2004 light curves, when the modulation 
index was 1.2% at 10.5 GHz and 0.5% at 4.9 GHz. 

The limited duration of the radio monitoring allows 
neither to fully constrain the true time scale of the 10.5 
GHz variations nor to check whether flux changes at 4.9 
GHz led or lagged those detected at the higher frequency. 
Hence, we are not able to say whether these flux changes 
are of intrinsic nature or if interstellar scattering is re- 
sponsible for at least part of the variations. Also a higher 
modulation index at higher frequency may be reconciled 
with interstellar scinti llation by invo king a strong scintil- 
lation regime (see e.g. iRickettl Il990|h 

However, in the following we speculate about the phys- 
ical implications of these variations in case we ascribed 
them to an intrinsic process. 

Under the assumption of an intrinsic nature of the 
10.5 GHz flux density variations detected in August 2, 
2004, the brightness temperature corresponding to the 
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Fig. 6. Radio light curves of the source and calibrators 
at 10.5 GHz (top, blue symbols) and 4.9 GHz (bottom, 
red symbols) obtained with the 100 m radio telescope at 
Effelsberg on January 28, 2005. The two calibrators are 3C 
067 (triangles) and 4C +09.11 (crosses). The horizontal 
lines indicate deviations of 1% from the mean flux den- 
sities, while the yellow areas highlight the period of the 
X-ray observations. 



initial fast decrease, calculated following IWagner et alJ 
l|l996|) 8 . is T b ~ 7 x 10 17 K, and a Doppler factor 6 > 
87 would be required to bring the rest frame value be- 
low the inverse Compton limit of 10 12 K (2~b )C .bserved = 
<5 3 Tb. intrinsic), confirming previous findings of a very 
high Doppler f actor to exp l ain some observing features 
of this source JKraus et all Ii999t iFuiisawa et all [l999; 
iFrev et allboOOHjorstad et alll200l[) . If we take the data 
errors into account, the inferred brightness temperature 
can actually vary from ~ 3 to ~ 11 x 10 17 K, correspond- 
ing to lower limits to S of 70 and 104, respectively. 

On January 28 the 10.5 GHz light curve shows a quasi- 
linear increasing trend, characterized by a ~ 2% flux vari- 
ation in about 4 hours. Under the assumption of intrinsic 
variability, the associated brightness temperature would 
be T b ~ 9 x 10 16 K, implying 5 > 45. Thi s value of the 
Doppl er factor is similar to that estimated bv lRaiteri et alJ 
(2005) from the Effelsberg observations of January 2004 
(6 > 46). Although it is much smaller than that derived 
from the August 2004 observations, nonetheless it would 
imply strong relativistic beaming of the emitted radiation 
in the intrinsic variability scenario. When considering the 

8 We adopt a fiat cosmology with Ho = 71 kms -1 Mpc -1 
and = 0.27, which yields a source luminosity distance of 
~ 6140 Mpc. 



data errors, the lower limit to 5 can actually range between 
23 and 62 for this variation. 

We have to mention that recent results, obtained dur- 
ing the current relatively quiescent state of the source, 
point to 8 and Tb values lower than the extrem e values es- 
timate d in the above mentioned earlier papers. IFrev et all 
l|2006|) analysed VSOP data acquired in 2001-2002. They 
found a decrease of the Doppler boosting by a factor of 
~ 20 with respect to the value S ~ 100 in ferred from 
data taken in 1999 just after a major outburst l|Frev et all 
20001). Their re sults are consistent with those derived by 
Kovalev et all l|2005|) from VLBA 2 -cm Survey data of 



March 20 01. Apparen t speed s up to (25.6 ± 7)c were mea- 
sured by IPiner et alJ l|2006h in the jet of AO 0235+164 
based on VLBA observations during 2002-2003. These 
val ues are quite smalle r than those previously reported 
bv Ijorstad et all l)200l[) and suggest Doppler factors not 
exceeding ~ 50. 

This boosting decrease can be explained, for instance, 
by an increased viewing angle of the rad io emitting region 
in the helical jet scenario proposed by IVillata fc Raiteri 
lll999l) and applied to the case of AO 0235+164 by 
lOstorero et al.l |2004). In this case, the high brightness 
temperature value we inferred from the August 2 data 
would be overestimated and the observed variation would 
be due, at least in part, to interstellar scattering. 

It is interesting to notice that the brightness tem- 
perature derived from VLBI observations performed on 
September 11, 2004 at 15 GHz is only 2.6 x 10 n K 
(MOJAVE collaboration, priv. comm.). Although both 
the VLBI and IDV brightness temperatures are lower lim- 
its, they are much more different th an it is usually foun d 
for these compact objects (see e.g. IZensus et all 12002^ . 
favouring an extrinsic explanation for the flux density vari- 
ation observed on August 2, 2004. 



5. Long-term behaviour 

The historical behaviour of the optical magnitudes and 
radio fluxes is plotted in Fig. These light curves have 
been shifted by a certain amount (indicated in the figure) 
to better distinguish the trend at different wavelengths. 
To increase the sampling in the optical, we plotted R- 
band data after JD = 2447000 and i?-band data converted 
into R ones before, according to R = B—<B — R>, with 
<B -R> = 1.7. 

The five yellow strips named a, b, c, d, and e high- 
light the radio outbursts (with optical count erpart) which 
were f ound to repeat quasi-periodically by iRaiteri et alJ 
ll200ll) . Indeed, the 1 yr thick strips, equally spaced by 
5.6 yr, contain the peak of the most prominent radio (and 
optical) outbursts. The last strip (f) indicates the period 
wh en the next ou t burst should have occurred according to 
the IRaiteri et al.l l)200l|) prediction, which motivated the 
WEBT campaign. 

It is evident that the last ~ 5 yr of suppressed activity 
in the radio bands are very atypical for this source. 
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Fig. 7. Historical optical (mag) and radio (37, 22, 14.5, 8, and 5 GHz, flux density in Jy) light curves of AO 0235+164 
shifted by the indicated amount; to increase the sampling in the optical, we show i?-band data after JD = 2447000 
and i?-converted S-band data before. Yellow strips highlight the periods (a, b, c, d, e, f) where a major radio (and 
optical) outburst is expected on the basis of a 5.6 ± 0.5 yr quasi-periodicity; periods a, (3, and 7 belong to an ~ 8 yr 
characteristic time scale of variability. 



If we run a periodicity analysis on the historical optical 
ligh t curves applyi ng the discrete Fourier transform (DFT; 
IPress et aD . [l9921. the most significant signal occurs at 
~ 3000 days (8.2 yr). A similar time scale can be found 
now also in the ra dio light curve s , toge ther with the 5-6 
yr period found bv lRaiteri et al.l (Jiboi) 9 . We show these 
results in Figs. E] and [5J 

Figure |H1 (top panel) displays the results of the DFT 
applied to the historical i?-band light curve. Before run- 
ning the DFT algorithm, data have been binned over 1 day 
intervals to homogenize the sampling. The significance of 
peaks in the DFT, which correspond to sinusoidal compo- 
nents of period r in the light curve, is estimated by the 
false alarm probability, i.e. the probability that a peak is 
higher than a certain level if the data are pure noise. The 
horizontal (dotted blue) line in the figure marks the value 
of power corresponding to a significance level better than 
0.001. The DFT indicates that the most likely period is 
8.2 yr, with other possible periods at 3.9 and 2.5 yr. The 
numerous signals corresponding to time scales around (or 
less than) 1 yr may be spurious, since they may reflect gaps 
in the data train. The bottom panel of Fig. El displa ys the 
autocorrelation funct i on (ACF: lEdelson fc Kroliki IlQsli 
iHufnagel fc Bregmarl Il992t IPeterson et all l2004h of the 
i?-band data. The most remarkable feature is that it shows 
two prominent peaks at ~ 7.6 and ~ 8.5 yr. The reason is 
that this method, which is more sensitive to maxima and 
minima brightness levels than to sinusoidal modulations 
of the light curves, feels the double-peaked nature of some 
events and the differences in the peak distances (see Fig. 
0). 

A similar plot for the 8 GHz light curve, the most 
extended and best sampled among the radio datasets, is 
presented in Fig. The DFT has been obtained after 
binning the data over 1 week, thus making the sampling 
more homogeneous. The strongest signal corresponds to 
the 5.6 yr period, which has not been suppressed by the 
lack of the predicted outburst in 2004, followed by signals 
at 3.7, 15.7, 8.0, 2.8, 1.9 yr, and by minor ones. We notice 
that some of the time scales are roughly multiples of the 
others, and that in the ACF plot the strongest peaks are 
found at about 11 and 15-17 yr, i.e. at twice the 5.6 and 



[years] 



We point out that the present analysis is performe d on 
different light curves with respect to lRaiteri et alJ (j2001): not 
only they are more extended in time, but also more datasets 
have been added to both the optical and radio light curves 
to improve the s ampli ng of the pre-2000 period, as discussed 
in iRaiteri et al] (12005^1 . In particular, besides the radio data 
from Metsahovi and UMRAO, we included data (both recent 
and past ones) from RATAN, Noto, Medicina, Effelsberg, and 
NRAO. 
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Fig. 8. Power spectrum of the discrete Fourier trans- 
form (top) and autocorrelation function (bottom) ob- 
tained from the historical i?-band light curve. In the top 
panel the blue dotted line indicates the level above which 
the significance is better than 0.001. 

8 yr periods detected with the DFT. However, even for 
these peaks the value of the ACF does not exceed 0.5. 

We also tested ou r resul ts with the structure function 
(SF; ISimonetti et"afl Il985h method, and found that its 
results are very similar to those already obtained with the 
ACF. 

An ~ 8 yr characteristic time scale of variability would 
put into correlation the outbursts occurred in the periods 
labelled a, /?, and 7 in Fig. [7| If this ~ 8 yr recurrence 
should repeat, an outburst will be visible during the 2006- 
2007 observing season. 

We also performed a cross-correlation analysis by 
means of the discrete correlation function (DCF) on the 
historical i?-band and radio light curves in order to check 
whether the addition of a further ob serving season of dat a 
can modify the results reported by IRaiteri et al. I ll2005lh 
Both the time lags of the radio variations with respect to 
the optical ones (from 1 to 2.5 months depending on the 
frequency separation) and the delays among variations at 
different radio frequen c ies (fr om 10 days to 1.5 months) 
found bv IRaiteri et al.l ll2005h are confirmed. 

6. Multiepoch spectral energy distribution 

The broad-band spectral energy distributions (SEDs) of 
AO 0235+164 corresponding to the times when X-ray 
satellites observed the source are displayed in Fig.^| Only 
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Fig. 9. Power spectrum of the discrete Fourier trans- 
form (top) and autocorrelation function (bottom) ob- 
tained from the historical 8 GHz light curve. In the top 
panel the blue dotted line indicates the level above which 
the significance is better than 0.001. 



epochs after 1990 were considered, thus neglecting the 
uncertain results obtained by Einstein in 1979-1981 and 
by EXOSAT in 1984. An analysis of the XMM-Newton 
and C handra observations can be found in IRaiteri et alJ 
(200(J, where a summary of previous X-ray observations 
by other satellites is also given together with the pertain- 
ing references. 

For the Chandra and XMM-Newton X-ray spectra, we 
display the results of both single power law (black bow- 
tie) and double power law (red line) models with Galactic 
absorption plus absor ption from the interv ening system at 
z — 0.524 obtained bv lRaiteri et alJ l)2006|) . Those authors 
favoured the latter, curved solution. For all epochs but the 
February 2002 one, for which the curved solution was by 
far statistically superior, this was mainly motivated by the 
fact that in this energy region we expect the superposition 
of the synchrotron and inverse-Compton components. We 
also display the UV and optical fluxes in the M2, Wl, 
U, B, and V bands derived by IRaiteri fit alJ feOOrfl by 
analysing the observations of the Optical Monitor onboard 
XMM-Newton. 

In the figure, optical (and near-IR) data simultaneous 
to the X-ray observations are plotted with filled symbols; 
when simultaneous data were missing, we searched for the 
closest data available, and showed them with empty sym- 
bols. Since radio variations are slower than the optical 
ones, we plotted as filled symbols both strictly simulta- 
neous radio data and data obtained as averages of those 
acquired just before and just after the satellite pointing. 

Optical data were corrected for the ELISA contribu- 
tion. Near-IR and optical magnitudes were de-reddened 
accord ing to the prescriptions given bv I Junkkajingn_et_al 



( 20041) . adopting the values tabulated bv IRaiteri et al 



(2005), and then tra nsformed into fluxes with the zero- 
point values given bv lBessell l|l979l) . 

Radio flux densities at 4.8, 8.0, and 14.5 GHz are 
from UMRAO, those at 22 and 37 GHz are from the 
Metsahovi Radio Observatory. RATAN simultaneous data 
at 0.9, 2.3, 3.9, 4.8, 7.7, 11.1, and 21.7 GHz contributed to 
some SEDs (ASCA-Feb98, Chandra- AugOO, XMM-Feb02, 
XMM-Aug04). In the last three SEDs, corresponding to 
the XMM-Newton pointings inside the WEBT campaign, 
radio data at 4.9 and 10.5 GHz from the Effelsberg tele- 
scope have been added. 

The green solid line plotted in each panel represents a 
cubic spline interpolation through the August 2004 data; 
it was drawn only to help the comparison of SEDs. 

As one can see from Fig.^J the flat radio spectrum can 
always be connected with the steep optical one through 
what is called "the synchrotron bump". This is likely 
peaking in the far-infrared, as low-energy peaked BL Lac 
objects are expected to do. Most of the times the X-ray 
spectrum is hard, suggesting that we are seeing the first 
part of the inverse-Compton bump, which is expected to 
peak in the GeV energy range. However, in summer 1993 
and February 2002 the X-ray spectrum is soft. As for the 
ROSAT spectrum, IMadeiski et all l)l996|) suggested that 
in a bright state the peak of the synchrotron component 
may shift to higher energies, and thus the synchrotron 
radiation would dominate in the soft X-ray region. Thi s 
interpretation was also followed bv lPadovani et al] l 2004h . 
when fitting the R OSAT spectrum with th e homogeneous, 
one-zone model of iGhisellini et alJ ll2002ft . However, this 
scenario is questioned by the contemporaneous optical 
data shown in Fig. Indeed, by looking at the vari- 
ous SEDs of this figure, there are no clues that the peak 
frequency may strongly change depending on the bright- 
ness state of the source. Actually, the situation may be 
more complex, since a UV-soft-X-ray bump, i.e. an ex- 
tra SED component in addition to the synchrotron and 
inverse-Compton ones, see ms to be present at lea st in 
some cases, a s discussed b y I.Tunkkarinen et all ((2004 and 
IRaiteri et all l|2005l EfJO^) . 

In the following, we briefly examine the individual 
SEDs, with particular attention to the optical data, since 
the presence of optical data simultaneous with the X-ray 
ones is fundamental to understand the emission mecha- 
nisms operating in the different energy bands. 

6.1. ROSAT: Summer 1993 

The ROSAT pointing lasted from July 21 to August 
15, 1993. The X-ray spectrum displayed as a blue, filled 
bo w-tie in the fig u re ref ers to the analysis performed 
by IMadeiski et al.l l|l996|) on the full dataset and re- 
veals a bright and soft X-ray spectrum. However, since 
the ROSAT light cu rve showed significant variability, 
IComastri et alJ 1 19971 reanalysed the data identifying two 
epochs of higher and lower count rates. The first epoch 
(July) yielded results similar to those of IMadeiski et alJ 
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Fig. 10. Spectral energy distributions of AO 0235+164 during X-ray satellite pointings. The green solid line represents 
a cubic spline interpolation through the August 2004 data. 



(1996). i?-band data taken on July 27 by the Hamburg 
Quasar Monitoring program 10 are plotted in the figure as 
blue triangles. The connection of these data with the X- 
ray spectrum seems to require an extra bump. 

The second ROSAT epoch (August) produced a dif- 
ferent spectrum, which is plotted in the figure as a black, 
empty bow-tie. It is simultaneous with the acquisition of 
the V and R data displayed as black dots, which were 
taken at the Lowell Observatory on August 13. In this 
case the extra bump is not strictly needed, but also not 
ruled out. 

6.2. ASCA: February 1994 and 1998 

The first ASCA pointing (bl ack, empty bow-tie) lasted 
from February 4 to 19, 1994 <|Madeiski et al.Ul996|) . The 
closest optical data (black circle s) were acqu i red a t the 
JKT (La Palma) on January 13 <|Takalo et al.Lll998() . 

ASCA pointed at the source again on February 
11-12, 1998, when it was in a bright radio and op- 
tical state (see Figs. 0] and Q). However, the X-ray 
flux was not particularly high, and the X-ray spectrum 
was harder (blue, filled bow-tie) than during the pre- 
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vious pointing ijjunkkarinen et all l2004|) . A simultane- 
ous optical-UV spectrum was acquired with the HST 
l|junkkarinen et ail 12004^ . which is displayed as a blue 
line in Fig. ^| and simultaneous BVR data were taken 
with the REOSC telescope of the Torino Observatory 
(blue squares), well matching the HST spectrum. The mis- 
alignment of the HST an d X-ray spectra was discussed by 
ljunkkarinen et all l|2004h . who considered the possibility 
of an extra SED component as unlikely and suggested in- 
stead that it may arise from the difficulties in modelling 
the extinction of the z — 0.524 absorber. 

It is interesting to point out here the correlation be- 
tween the radio and X-ray spectral states: the harder and 
brighter X-ray spectrum of 1998 corresponds to a brighter 
and flatter radio spectrum, as expected when considering 
the inverse-Compton origin of the X-ray radiation. 

6.3. RXTE: November 1997 



The RXTE pointing was performed on November 3-6, 
1997, during the rising phase of the o ptical outburst (se e 
Fig. EI). The X-ray flux was found low l|Webb et alll2000|) . 
Simultaneous optical data were taken by many observa- 
tories in the ambit of the first-light WEBT campaign 
[ Raiteri et all l200lj) . In Fig. QUI we plotted two BVRI 
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spectra obtained at the Teide Observatory on November 
3 and 6 in order to have homogeneous data. Also shown 
is a point at 850 pm obtained by averaging data taken 
at the JCMT befo re and after the RXTE observations 
llWebb et alll200rt) . 

6.4. BeppoSAX: January 1999 

BeppoSAX observed AO 0235+164 on January 28, 1999. 
The bow-tie in Fig. II 0| shows the results of the analysis by 
Padovani et alJ |2004), while the points in the X-ray do- 
main were obtained directly by the ASI BeppoSAX Data 
Center 11 . The i?-band datum was taken at the Torino 
Observatory 5 days before. 

6.5. Chandra: August 2000 

The Chandra observations were performed on August 20- 
21, 2000. They reveale d a faint and hard X -ray spectrum 

I Turnshek et~afl 120031: iRaiteri et all kOOd) . The closest 
optical datum was taken at the Torino Observatory 7 
weeks later. 

6.6. XMM-Newton: February 2002 

The XMM-Newton o bservations of Febru ary 10, 20 02 were 
analysed in det ail by Raite ri et al ] l|200flh (see also lKadlerl 
2005] iFoschini et alJl200fih . The source was found in the 
brightest X-ray state ever detected (Fx koV ~ 2 p Jy) , with 
a soft spectrum, and the X-ray light curve revealed no- 
ticeable variability on a time scale of ~ 40 minutes. The 
Optical Monitor (OM) onboard XMM-Newton observed 
only in the UVIU1 band, finding a bright UV state. Three 
i?-band frames were acquired at the Mount Maidanak 
Observatory ~ 1.5 days before the XMM-Newton obser- 
vations, revealing a moderately bright state (R ~ 16.6). 
In the figure we plotted a complete BVRI optical spec- 
trum taken at Mount Maidanak 11 days earlier, when the 
source was exactly at the same optical level. However, an 
/-band datum taken at the Torino Observatory 9 days 
after the XMM-Newton pointing showed the source to 
be much fainter (/ ~ 16.9), so it is possible that at the 
time of the XMM-Newton observation the source was in a 
dimming phase, and that simultaneous optical data would 
place just below those shown in the figure. In conclusion, 
it is very unlikely that the optical, UV, and X-ray data can 
be fitted by a single convex synchrotron component. The 
UV-soft-X-ray bump would instead provide a comfortable 
way to describe this SED. 

6. 7. XMM-Newton: January 2004 

The SED corresponding to the XMM-Newton point- 
ing of Janua r y 18- 19, 2004 was described in detail in 
IRaiteri et alJ l)2005|) . The contemporaneous radio data 
come from Metsahovi, Effelsberg, and UMRAO, the op- 

II http : //www. asdc . asi . it/bepposax/ 



tical (and near-IR) data were taken at the NOT, while 
the OM observed in the V, B, U, UW1, and UVA/2 fil- 
ters. The transition from a steep optical spectrum to a flat 
UV one is dramatic, and the probable curvature of the X- 
ray spectrum supports the existence of the UV-soft-X-ray 
bump. 

6.8. XMM-Newton: August 2004 

The plotted radio points at 4.9 and 10.5 GHz are aver- 
ages of the intraday observations at Effelsberg (see Fig. 
El; the 14.5 GHz datum is from UMRAO; the 37 GHz 
point represents the average of the Metsahovi data. Data 
at 2.3, 4.8, 7.7, 11.1, and 21.7 GHz taken simultaneously 
by RATAN at the beginning of September have been 
added to increase the sampling. They nicely agree with the 
other radio data. Near-IR information comes from Campo 
Imperatore. Many optical data are available in the BVRI 
bands, with a fair sampling, showing a constant level with- 
out any significant variation. We display the averages of 
the data taken at Skinakas, to have a homogeneous BVRI 
set. Other BVRI sequences acquired at Mount Maidanak 
just before the start of the XMM-Newton observations 
confirm the Skinakas average spectrum. The hints for the 
extra bump come from the UVA/2 datum and the strong 
curvature of the X-ray spectrum when a double power law 
model is applied. 

6.9. XMM-Newton: January 2005 

As in the previous epoch, the points at 4.9 and 10.5 GHz 
represent the averages of the Effelsberg data (see Fig.EJ; 
the 14.5 GHz datum is from UMRAO, the 37 GHz one 
was obtained at Metsahovi, and the 2 and 3 mm data 
were acquired at Pico Veleta. In the optical band, the 
/ datum is from Mount Maidanak and the R datum is 
from Boltwood. Here the curvature of the X-ray spectrum 
is smaller than in the previous epoch, but there are two 
UV points (UVW1 and UVAf 2 from the OM) that may 
support the existence of a UV bump. 

In summary, the existence of the UV-soft-X-ray bump 
in the SED of AO 0235+164 is supported by two argu- 
ments: the SED hardening in the UV domain, which gives 
rise to a concave optical-UV branch, and the correspond- 
ing "rise" of the soft-X-ray spectra towards the UV band. 
Even if the first argument should fall in case the UV ex- 
tinction in the z = 0.524 intervening system had been 
strongly overestimated bv ljunkkarinen et alJ l|2004|) . the 
extra-bump cannot be avoided in the cases where the X- 
ray spectrum is soft. 

7. Discussion and conclusions 

Among blazars, AO 0235+164 presents several peculiar 
features. It has shown dramatic flux changes at all wave- 
lengths on a wide range of timescales, including IDV. Its 
main outbursts have repeated quasi-periodically, at least 
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for some time, with clear correlation of the radio and op- 
tical flux variations with very short time delays. When 
considering this along with the presence of an intervening 
galactic system, also gravitational microlensing becomes a 
viable interpretation for the source variability. At milliarc- 
second scales it shows an extremely compact core, with 
faint extensions occasionally appearing at various posi- 
tion angles. Several observations suggest that the source 
radiation is unusually highly Doppler-boosted. 

In this paper we presented the results of the radio, 
near-IR, and optical monitoring during the second observ- 
ing season of the WEBT campaign 2003-2005 targeting 
AO 0235+164. The results of the first observing season 
and a detailed analysis of the three XMM-Newton obser- 
vations per formed during the campaig n have already been 
reported in lRaiteri et all l)2005t |2006|) . 

The big radio outburst, with a possible optical counter- 
part, which was foreseen to peak around February-March 
2004 on the basis of the quasi-p eriodicity of 5.7 ± 0.5 yr 
discussed bv lRaiteri et al.ll|200l|h did not occur. The radio 
flux remained rather low during the whole WEBT cam- 
paign, with some mild increased activity at the higher 
radio frequencies only in the last observing season 2004- 
2005. A mild increased activity in the last season also 
characterized the optical and near-IR light curves. The 
comparison of the optical and radio data taken during the 
WEBT campaign with the long-term light curves of this 
object shows how this last, long period of suppressed ac- 
tivity especially at radio wavelengths is unusual for this 
source. The analysis of the historical light curves reveals 
now a longer characteristic time scale of variability of ~ 8 
yr, more evident in the optical, but also present in the 
radio data along with the ~ 5.7 yr one. 

The colour-index analysis reveals noticeable spectral 
changes in the optical band, which nonetheless do not 
show correlation with the source magnitude. The long- 
term trend appears essentially achromatic, with a possible 
progressive reddening in the last years. 

All the three XMM-Newton observations performed 
during the WEBT campaign found the source in a faint 
state at all wavelengths. Neither the optical nor the X-ray 
light curves show any significant variability. In contrast, 
the simultaneous intraday radio observations at the 100 m 
Effclsberg telescope reveal some flux changes at 10.5 GHz. 
If they are ascribed to intrinsic processes, a Doppler factor 
up to ~ 100 can in principle be derived, confirming pre- 
vious findings that this is an extremely beamed source. 
However, on the basis of these data we cannot rule out 
that interstellar scintillation was responsible for at least 
part of the variability. Actually, a reduced Doppler fac- 
tor would agree with recent findings pointing to a reduced 
beaming during the current relatively quiescent period. 

We compiled the SEDs corresponding to the X-ray 
observations since the ROSAT pointing in summer 1993. 
Radio and optical data simultaneous with the X-ray ob- 
servations have been looked for in our database in or- 
der to shed light on the various emission components in- 
volved. The aim was mainly to find clues to the pres- 



ence of the UV-soft- X-r ay b ump already discussed b y 
lJunkkarinen et alJ l|2004[) and lRaiteri et all (|2005l 12006^ . 
The existence of this extra component, besides the classi- 
cal synchrotron and inverse-Compton ones, is mainly sup- 
ported by the UV data, when available, and by the soft 
X-ray spectra detected by ROSAT in summer 1993 and 
by XMM-Newton in February 2002. 

A discussion on th e possible int e rpreta tions of this 
bump can be found in iRaiteri et al.l (|2006). Among the 
various pictures, these authors focused mainly on the pres- 
ence of a thermal accretion disc and on an additional 
synchrotron component. In the latter case, which was 
favoured by the authors, we can envisage two possibil- 
ities. Either there are actually two synchrotron compo- 
nents, coming from distinct and separate jet regions with- 
out any simple relationship, or we are in the presence of a 
continuous inhomogeneous jet whose radiation is for some 
reason suppressed in the optical-near-UV region. Here we 
have again two possibilities: i) the break in the spectrum 
is due to a discontinuity in the o pacity (as in the case of 
Mkn 501. lRaiteri fc Villatal20o"3j) . ii) the suppressed emis- 
sion is caused by a misalignment of the corresponding jet 
region whose radiation is less beamed towards the line of 
sight 12 . The observed strong variability in both these SED 
components might favour a scenario where both the emis- 
sion from the near-IR to the soft X-rays and the spectral 
gap change in time as the strongly bent jet precesses or 
rotates. This view may be supported by the large vari- 
ations of the VLBI jet compone nts position angles ob- 
served by many authors (see e.g. iFrev 
l.lorstad et all l200li IPiner et allfeoorf) . 

We noticed how among the ten satellite pointings in 
1993-2005, only two of them (ROSAT in summer 1993 
and XMM-Newton in February 2002) found the source in a 
bright and soft X-ray state, possibly of synchrotron origin. 
The All Sky Monitor X-ray light curve of AO 0235+ 164 13 
from January 1996 to November 2005 shows that X-ray 
bright states are indeed very uncommon for this source: 
only in the first half of 2002 a significant increase of the 
count rate was observed, in correspondence to the first 
XMM-Newton pointing. 

The two bright X-ray states are separated by about 
8.5 yr, and a ~ 9 yr interval elapsed between the ROSAT 
and the 1984 EXOSAT observation, when a high (but very 
uncertain) X-ray flux was measured. A similar time sepa- 
ration is found when considering the last strongest optical 
outbursts of 1990 and 1998. Thus, both of them seem to 
follow an X-ray brightening by about 5 yr. It will be in- 
teresting to see whether an optical flux increase will occur 
in the next observing seasons, i.e. about 5 yr after the 
2002 X-ray bright state. It would support both the possi- 
ble quasi-periodicity of ~ 8 yr of the optical (and radio) 



12 Although at other wavelengths (and hence sp atial resolu- 
tion) , this scenario has indeed been suggested by iBach et alJ 
(2006) to interpret the region of suppressed flux density in the 
intensity profiles of the VLBA images of BL Lacertae. 

13 http : //heasarc .nasa. gov/docs/xte/ASM/ sources .html 
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outbursts and a correlation between the optical and X- 
ray emission with a long time delay. A correlation of this 
type would be explained by a scenario where a curved in- 
homogeneous jet precesses or rotates, thus subsequently 
aligning its different emitting portions with the line of 
sight. 

Several WEBT members are continuing to monitor the 
source in the radio, near-IR, and optical bands so that we 
will be able to check whether the source activity continues 
to increase and eventually leads to a forthcoming outburst. 
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